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Abstract 

We compute the branching ratios for fj, ej, t ^ fi^ and r — > 67 in a super- 
symmetric model invariant under the flavour symmetry group A4 x Z3 x U{1)fn, in 
which near tri-bimaximal lepton mixing is naturally predicted. At leading order in 
the small symmetry breaking parameter u, which is of the same order as the reactor 
mixing angle 613, we find that the branching ratios generically scale as u^. Applying 
the current bound on the branching ratio of /i — > 67 shows that small values of u or 
tan/3 are preferred in the model for mass parameters msuSY and mi/2 smaller than 
1000 GeV. The bound expected from the on-going MEG experiment will provide a 
severe constraint on the parameter space of the model either enforcing u ~ 0.01 and 
small tan/3 or msusY and mi/2 above 1000 GeV. In the special case of universal soft 
supersymmetry breaking terms in the flavon sector a cancellation takes place in the 
amplitudes and the branching ratios scale as tt^, allowing for smaller slepton masses. 
The branching ratios for r — > and r — > 67 are predicted to be of the same order 
as the one for fi 67, which precludes the possibility of observing these r decays in 
the near future. 
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1 Introduction 



Flavour violation in the lepton sector (LFV) has been firmly established in neutrino oscil- 
lations. A direct evidence of conversion of electron neutrinos into a combination of muon 
and tau neutrinos is provided by solar neutrino oscillations. In a three neutrino framework 
atmospheric neutrino data can only be explained by the dominant oscillation of muon neu- 
trinos into tau neutrinos. It is natural to expect that LFV takes place, at least at some 
level, in other processes such as those involving charged leptons. Flavour violating decays 
of charged leptons, strictly forbidden in the Standard Model (SM), are indeed allowed as 
soon as neutrino mass terms are considered. If neutrino masses are the only source of 
LFV, the effects are too small to be detected, but in most extensions of the SM in which 
new particles and new interactions with a characteristic scale M are included, the pres- 
ence of new sources of flavour violation, in particular in the lepton sector, is a generic 
feature. In a low-energy description, the corresponding effects can be parametrized by 
higher- dimensional operators describing flavour-violating rare decays of the charged lep- 
tons. The dominant terms are represented by dimension six operators, suppressed by two 
powers of M. The present bounds on the branching ratios [1] of these decays set strin- 
gent limits on combinations of the scale M and the coefficients of the involved operators. 
Typically, for coefficients of order one, the existing bounds require a large scale M, several 
orders of magnitude larger than the TeV scale. Conversely, to allow for new physics at 
the TeV scale, coefficients much smaller than one are required, which might suggest the 
presence of a flavour symmetry. 

New physics at the TeV scale supplemented by a flavour symmetry provides an in- 
teresting framework. Flavour symmetries have been invoked to describe the observed 
pattern of lepton masses and mixing angles [2], but quite often this approach is limited 
to a flt of the existing data, with very few new testable predictions [3]. Speciflc relations 
among LFV processes are usually consequences of flavour symmetries and of their pat- 
tern of symmetry breaking. This allows to get independent information on the flavour 
symmetry in the charged lepton sector. Moreover, if M is sufficiently small, new parti- 
cles might be produced and detected at the LHC, with features that could additionally 
confirm or reject the assumed symmetry pattern. In [4] we have recently analyzed in an 
effective Lagrangian approach LFV processes within a class of models possessing a flavour 
symmetry Gf = A4 x x U{1)fn- Models of this class [5-8] automatically reproduce 
nearly tri-bimaximal (TB) lepton mixing at the leading order (LO) via a vacuum alignment 
mechanism. In this type of models [5-8, 10, 11], corrections to TB mixing are generically 
proportional to the symmetry breaking parameter u, especially the reactor mixing angle 
6'i3 is of the order of The parameter u is expected to vary between a few per mil and 
a few percent. In particular we have evaluated the normalized branching ratios Rij for the 
LFV transitions k — > Ij'y: 

BR{k-.l,^) 
BR{k^ljViVj) ■ ^' 

^Due to the success of these models in the lepton sector several extensions to the quark sector can be 
found in the literature [9]. 

^For alternative scenarios see [12,13]. 
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In [4] we found that the generic expectation for the ratios Rij is: 

= ^^p^Fii «l (2) 

where a is the fine structure constant, is the Fermi constant and Wij are dimensionless 
parameters that cannot be predicted within the effective Lagrangian approach, but are 
expected to be of order one. Given the range of m, it turns out that the present bound on 
/i — > 67 requires the scale M to be larger than about 10 TeV. If the underlying fundamental 
theory is weakly interacting, this translates into a lower bound of about gM/ (47r) on the 
typical mass of the new particles. At variance with this generic estimate, it was also 
observed that, if the flavour symmetry is realized in a supersymmetric (SUSY) model, a 
cancellation might take place in the amplitudes for LFV transitions. An argument was 
given suggesting that in this case the ratios RfJ^^"^ are of the form: 



SUSY ^ 487r^a 



I (1) 2|2 , (2) I 



(3) 



where M now corresponds to the supersymmetry breaking scale, rrii {i = e, /i, r) are the 
charged lepton masses and w^j ' are unknown dimensionless quantities of order one. There- 
fore in a supersymmetric context the predicted rates are more suppressed, allowing for new 
physics close to the TeV scale, without conflicting with the present bounds. Since the can- 
cellation expected in the supersymmetric case considerably changes the conclusion that can 
be derived from the existing bound on /i — >■ 67, it is important to perform a direct compu- 
tation of the branching ratios within an explicit supersymmetric model incorporating the 
flavour symmetry x Z3 x U{1)fn- 

In the present paper we consider the A4 realization proposed in [4] in which SUSY 
breaking effects were ignored and we extend it to a more realistic model, by adding a full 
set of SUSY breaking terms consistent with the flavour symmetry. We assume that the 
breaking of supersymmetry occurs at a scale higher than or comparable to the flavour scale, 
simulated in our effective Lagrangian by a cutoff, so that at energies close to the cutoff scale 
we have non-universal boundary conditions for the soft SUSY breaking terms, dictated by 
the flavour symmetry. Depending on the assumed mechanism of SUSY breaking we may 
have boundary conditions different from these, possibly enforced at a smaller energy scale. 
For this reason, our approach maximizes the possible effects on LFV processes. 

Through a detailed calculation of the slepton mass matrices in the physical basis, in 
which kinetic terms are canonical and leptons are in their mass eigenbasis, and evaluation 
of Rij in the mass insertion (MI) approximation we find that the behaviour of Rij, expected 
from the SUSY variant of the effective Lagrangian approach, given in eq. ^ is violated by 
a single, flavour independent contribution. The correct scaling of Rij is the one of eq. (|2]), 
with a universal constant Wij. This implies i?^e = Rrfi = Rre at the LO in u. We identify 
the source of violation of the expected behaviour in a contribution to the right-left (RL) 
block of the slepton mass matrix, associated to the sector necessary to maintain the correct 
breaking of the flavour symmetry A4. We also enumerate the conditions under which such 
a contribution is absent and the behavior in eq. ^ is recovered, though we could not find 
a dynamical explanation to justify the realization of these conditions in our model. 
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As shown in the MI approximation the coefficient Wij of the leading contribution in eq. 
([2]) is numerically small in a large region of the parameter space, where the sub-leading 
contributions play an important role. We then provide a numerical study of Rij by using full 
one-loop expressions and explore the parameter space of the model. Thereby, we assume 
a supergravity (SUGRA) framework with a common mass scale msusY for soft sfermion 
and Higgs masses and a common mass mi/2 for gauginos at high energies. Applying 
the current MEGA bound on i?^e we find that small values of u or tan/3 are favoured 
for msusY and mi/2 below 1000 GeV. Employing the foreseen bound coming from the 
MEG experiment constrains the model severely allowing only for small u and tan/5 for 
msusY, ^1/2 ^ 1000 GeV or requires larger msusv and mi/2 reducing the prospects for 
detection of sparticles at LHC. Furthermore, it turns out to be rather unnatural to reconcile 
the values of superparticle masses necessary to account for the measured deviation 6a^ in 
a^, the muon anomalous magnetic moment, from the SM value with the present bound on 
Rfj^e- In our model values of (5a^ smaller than 100 x 10~^^ are favoured. 

The paper is structured as follows: in section 2 we review the basic features of the 
SUSY A4 model and calculate kinetic terms and mass matrices of leptons and sleptons in 
an expansion in the symmetry breaking parameters. In section 3 we compute the slepton 
mass matrices in the physical basis. Furthermore, we estimate the possible effects on the 
slepton masses coming from renormalization group (RG) running. Section 4 is dedicated 
to the study of the quantities Rij in the MI approximation, while in section 5 we discuss 
Rfj,e in a more quantitative way in the SUGRA context and comment on the size of the 
measured deviation in the muon anomalous magnetic moment in our model. Finally, 
we conclude in section 6. Appendix A contains details about the group theory of A4, 
appendix B details of the calculation of the slepton mass matrices in the physical basis 
and in appendix C conventions and formulae, used in the numerical study of section 5, are 
given. 

2 A SUSY model with A4 flavour symmetry 

In this section, we ffist discuss the main features and predictions of the lepton mixing in 
the class of ^4 realizations to which the present SUSY model belongs. Then we will present 
the Kahler and the superpotential of our model from which kinetic and mass terms are 
derived for leptons and sleptons. 

2.1 General features 

The flavour symmetry A4 can give rise to TB mixing in the lepton sector, if it is broken in 
a specific way, as discussed in [5-8, 10]. The group A4 is the group of even permutations 
of four objects which is isomorphic to the tetrahedral group. It has 12 elements and four 
irreducible representations: three inequivalent singlets, denoted by 1, 1' and 1" and a triplet 
3. A4 is generated by two generators S and T which fulfill the following relations [14]: 

52 = {STf = T^ = 1 . (4) 
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A set of generators S and T for all representations as well as Kronecker products and 
Clebsch Gordan coefficients are given in appendix A. The specific breaking pattern of A4 
which leads to TB mixing for leptons requires that a Z3 subgroup, called Gt, is preserved 
in the charged lepton sector, whereas a Z2 group, Gs, remains unbroken in the neutrino 
sector. The Z3 group is generated by the element T and is left intact, if a scalar triplet 
ip = (p2, fs) acquires a vacuum expectation value (VEV) of the form: 

(<^)oc (1,0,0). (5) 

Breaking ^4 down to 6*5 is achieved by a triplet which gets a VEV: 

(^)oc (1,1,1). (6) 

Apart from A4, being responsible for TB mixing, the theory is invariant under additional 
fiavour symmetries so that the full symmetry is: 

G'^ = A4xZ3xU(l)^^xU(l)^. (7) 

The cyclic symmetry Z3 (not to be confused with the Z3 subgroup of ^4 preserved at LO 
in the charged lepton sector) is necessary in order to avoid large mixing effects between the 
fiavons that give masses to the charged leptons and those giving masses to neutrinos. The 
mass hierarchy among the charged leptons is explained through the U(l)^^ factor [15]. The 
second U(l) factor, U(l)^, is a continuous i?-symmetry containing i?-parity as a subgroup. 
The flavour symmetry Gf is broken by the chiral superflelds ipT, ^ and 6fn whose 
transformation properties under Gf are shown in table 1 together with those of the lepton 
supermultiplets /, e'^, ^f, and of the electroweak doublets Hu^. Introducing the driving 
flelds v?o 5 V^o ^0, see also table 1, we can write down the following superpotential: 

Wd = Mt{(pI(Pt) + gi^l^T^x) 

+ 9i{vlvsVs) + 92ii^o^s) + Qs^oi^s^s) + Qiioi'^ + Qbioil + g&iol'^ + •••• (8) 

where dots denote sub-leading non-renormalizable corrections. By (■ ■ • ) we denote the 
contraction to an invariant. In the limit of unbroken supersymmetry all F terms 
vanish. From these conditions we can derive the vacuum alignment of the flavons (/9t, '^s 
and ^: 



A/ 



(m, 0, 0) + (cV, CM^ cu^) + 0{u^) 

Ck{u,u,u) + 0{u^) (9) 

CaU + 0{U^) 



A/ 
A/ 

where c, c', Ca^b are complex numbers with absolute value of order one and u is one of the 
two small symmetry breaking parameters in the theory. The scale Aj is the cutoff scale 



■^Thc additional flavon field ^ transforms in the same way under the symmetries of the model as the 
field 1^. However, it does not acquire a VEV at LO and thus is not relevant for our discussion here. For 
further details see [6]. 
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associated to the flavour symmetry. Its value is expected to be around the scale of grand 
unification 10^^ GeV. We display the sub-leading O(m^) corrections, which arise from non- 
renormalizable terms in the superpotential and which have been studied in detail in [6] , 
only for the field ifT, since only these are relevant in the following analysis, whereas those 
of the fields ^ and ips are important in the computation of neutrino masses and lepton 
mixings. We remark that in the SUSY limit the VEVs of the driving fields ip^, if^ and 
are zero. This is however in general no longer true, if we include soft SUSY breaking terms 
into the fiavon potential, as we shall see in the next section and as has been discussed 
in [16]. 

For the Froggatt-Nielsen (FN) field dpNto acquire a VEV, we assume that the symmetry 
U(l)^^ is gauged such that 9fn gets its VEV through a D term. The corresponding 
potential is of the form: 

hM^^j - gp^iep^]-^ + ...y (10) 



Vd. 



FN 



where gpN is the gauge coupling constant of U(l)p^ and Mpj denotes the contribution 



of the Fayet-Iliopoulos (FI) term. Dots in eq.(lO) represent e.g. terms involving the 
right-handed charged leptons e'^ and /i"^ which are charged under U(l)p,^. These terms are 
relevant in the calculation of the different contributions to the soft mass terms of right-right 
(RR) type, as we shall see below. Vd,fn leads in the SUSY limit to: 



\{0 



FN) 



which we parametrize as: 



'FAT 



) 



A 



9fn 



t 



f 



(11) 



(12) 



with t being the second small symmetry breaking parameter in our model. Both u and t 
can also be in general complex, but through field redefinitions of the supermultiplets ipx 
and 9fn, they can be made real and positive. 



Field 










Hu,d 










T 




^0 


A, 


3 


1 


I" 


V 


1 


3 


3 


1 


1 


3 


3 


1 


Z3 


UJ 








1 


1 


UJ 


UJ 


1 


1 


UJ 


UJ 


U(l)p^ 





2 


1 

















-1 











U(l)fi 


1 


1 


1 


1 

















2 


2 


2 



Table 1: Transformation properties of the fields under A^, Z3, U(l)^^ and U(l)^. 



At the LO, the mass matrices of charged leptons and light neutrinos have the following 
form: 

Vet^ 

mi(x \ y^t I M , (13) 

Vr 
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and 

/ a + 26/3 -6/3 -6/3 \ 

oc -6/3 26/3 a - 6/3 u , (14) 

\ -6/3 a - 6/3 26/3 / 

with ?/e, y^, Ur, a and 6 being complex numbers with absolute values of order one. As one 
can see the relative hierarchy among the charged lepton masses is given by the parameter 
t for which we take: 

t^Om . (15) 
The neutrino mass matrix is diagonalized by: 

U^^m^UTB oc diag(a + 6, a, —a + b)u , (16) 

where Utb, up to phases, is the TB mixing matrix [17]: 

/ ^/2/3 1/V3 \ 
Utb = -1/V6 1/V3 -1/v^ . (17) 
V -1/V6 1/V3 +1/V2 J 

The allowed range of the parameter u is determined by the requirement that sub-leading 
corrections which perturb the LO result are not too large and by the requirement that the 
T Yukawa coupling does not become too large. The first requirement results in an upper 
bound on u of about 0.05, which mainly comes from the fact that the solar mixing angle 
should remain in its la range [2]. The second one gives a lower bound which we estimate 
as: 

tciia(3y/2mr tan/5 
u = - — . ^ 0.01- — - (18) 

where v ~ 246 GeV and tan/5 is the ratio between the VEVs of the neutral spin zero 
components of Hu and Hd, the two doublets responsible for electroweak symmetry breaking. 
For the r lepton we use its pole mass = (1776.84 ±0.17) MeV [18]. Requesting \yr\ < 3 
we find a lower limit on u close to the upper bound 0.05 for tan /3 = 15, whereas tan f3 = 2 
1^ gives as lower limit u > 0.007. Obviously, these limits depend on the largest allowed 
value of \yr\, as well as on whether we identify with the pole mass or with the r mass 
renormalized at some reference scale, such as the scale of grand unification. We choose as 
maximal range: 

0.007 <n<0.05 , (19) 

which shrinks when tan/5 is increased from 2 to 15. Concerning the relative size of the two 
symmetry breaking parameters we note that m < t holds for all values of u. 

Ht is known that tan /? cannot be too small [19]. Here we take tan f] — 2 as the smallest allowed value. 
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2.2 The SUSY Lagrangian and the soft SUSY breaking terms 



We analyse the Lagrangian of the modeh 

C = j (f9(fe}Ciz,e^^z)+ j d^Owiz) + h.c. 

(fef{z)WW + h.c. 



+ 



(20) 



where }C{z, z) is the Kahler potential, a real gauge-invariant function of the chiral super- 
fields z and their conjugates, of dimensionality (mass)^; w{z) is the superpotential, an 
analytic gauge-invariant function of the chiral superfields, of dimensionality (mass)^; f{z) 
is the gauge kinetic function, a dimensionless analytic gauge-invariant function; V is the 
Lie-algebra valued vector supermultiplet, describing the gauge fields and their superpart- 
ners. Finally W is the chiral superfield describing, together with the function f{z), the 
kinetic terms of gauge bosons and their superpartners. Each of the terms on the right- 
hand side can be written in an expansion in powers of the flavon fields. Since we have 
two independent symmetry breaking parameters, u and t, see eqs. Q and eq. (Q, we 



consider a double expansion of C in powers of u and t. In this expansion we keep terms up 
to the second order in u, i.e. terms quadratic in the fields ips,T and ^. The expansion in 
the parameter t, responsible for the breaking of the Froggatt-Nielsen l]{l)pj\f symmetry, is 
stopped at the first non-trivial order, that is by allowing as many powers of the field 9fn 
necessary in order to obtain non-vanishing values for all entries of the matrices describing 
lepton masses as well as for the entries of the matrices describing kinetic terms and slepton 
masses. Finally, second order corrections in u also arise from the sub-leading terms of 
the VEV {^t), eq. Q, and are included in our estimates. 

The soft SUSY breaking terms are generated from the SUSY Lagrangian by promoting 
all coupling constants, such as Yukawa couplings, couplings in the flavon superpotential and 
couplings in the Kahler potential, to superfields with constant 6^ and 6^6 components [20] . 
Through this we derive subsequently the soft masses ijn\ev)LL)K a^iid ijnlp^j^K from the 
Kahler potential. One contribution to Tnl^^^, which we call (^eijL)i following, arises 

from the Yukawa couplings present in the superpotential w. 

Important contributions to slepton masses originate from the modification of the VEVs 
of flavons and driving fields due to SUSY breaking effects. A detailed study of the VEVs 
of these fields and their dependence on the soft SUSY breaking parameters is presented 
in [16] and we summarize the main results here. When soft SUSY breaking terms are 
included into the flavon potential, the VEVs in eq. ^ receive additional contributions 
of order msusY, completely negligible compared to A/ m. At the same time, the driving 
fields (fQ, and develop a VEV of the size of the soft SUSY breaking scale msusY- 
An equivalent statement is that the auxiliary fields of the fiavons acquire a VEV at the 



^Concerning the Kahler potential we observe that we can additionally write down operators involving 
the total invariant 0fn0fn = \0fn\'^ ■ These contribute to the diagonal elements of the kinetic terms and 
the slepton masses. In the Kahler potential for the left-handed fields they can be safely neglected, since 
the LO correction is of 0{u). In the right-handed sector, they contribute at the same order as the terms 
arising through a double flavon insertion. 
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LO of the size of msusY x uAf. Especially, for the auxiliary fields contained in the fiavon 
supermultiplet ipr we have [16]: 

C msusY { {u, 0, 0) + {c'pu^, cpu^, cpu^) } (21) 



A J \ d(f 

where dp and cp are in general complex numbers with absolute value of order one. 
The parameter C, vanishes in the special case of universal soft mass terms in the fiavon 
potential. When different from zero, the VEVs of the auxiliary components of the fiavon 
supermultiplet (pr generate another contribution to the soft masses of RL-type, which 
we denote as {jn1^p)2. This contribution is analogous to the one which has been found 
before in the supergravity context and which can have a considerable effect on the size of 
the branching ratio of radiative leptonic decays, as shown in [21,22]. Indeed, as we shall 
see below, in the global SUSY model under consideration the leading dependence of the 
normalized branching ratios Rij on u is dominated by {fn^Rhl'i- We remark that the VEVs 
in eq. (21 ) and those of the corresponding fiavon field i^t in eq. ^ have a similar structure 



but they are not proportional, in general. This is due to the different coefficients c, c' and 
Ci?, dp, which can be qualitatively understood as follows: the coefficients c, d mainly 
depend on a set of parameters that remain in the SUSY limit and receive completely 
negligible corrections from the SUSY breaking terms. On the contrary {dw/dipT) vanishes 
in the SUSY limit, to all orders in m, and cp, dp crucially depend on the set of parameters 
describing the SUSY breaking. We will see that, if c and cp accidentally coincide (up to 
complex conjugation), a cancellation in the leading behaviour of Rij takes place. 

Similarly, the VEV of the FN field 6fn becomes shifted, when soft SUSY breaking 
terms are included into the potential, so that: 

/l/f2 

-^-\{epN)? = ceml^SY^ (22) 
9fn 

with C0 being an order one number, holds. This will lead to a contribution {jn1pi^j^£,,FN 
to the soft masses of RR-type, since only the right-handed charged leptons d" and ijf^ 
are charged under U(l)p^. Apart from these there are supersymmetric contributions 
to Tn^ej^-jpp and mg^^ from F and D terms, {rn'^eu)LL)F{D) and ijnl^^pi^p,), as well as a 
contribution to ml^p coming from the F term of Hd-, called {ml^i)^ in the following. 

In our notation a chiral superfield and its -R-parity even component are denoted by the 
same letter. The i?-parity odd component is indicated by a tilde in the following and the 
conjugate (anti-chiral) superfield is denoted by a bar. 



2.2.1 Kahler potential 

The expansion of the Kahler potential can be written as: 

/C = /C(°) + /C« + /C(2) + ... (23) 

In our model the Kahler potential deviates from the canonical form, )C{z,z) = zz, due 
to the contributions of non-renormalizable terms, invariant under both the gauge and the 
flavour symmetries. Such contributions are sub-leading in the {u, t) expansion, but in 
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principle they are non-negligible, since the redefinitions required to arrive at canonically 
normalized fields can also affect lepton and slepton mass matrices. The LO term in u is 
given by: 



l<'^l^ + K\HX + kd\H£ + kFN\OFN\^ + .■. (24) 



where l'^ = (e^, /i'^, t'^) and dots stand for additional contributions related to the flavon sec- 
tor. The quantities kQ, (/cq),, (A;o)j, k^^d, kpN and the corresponding ones in the following 
formulae, are treated as superfields with a constant O'^O component, as is explained above. 
Concerning the effect of the superfield OpN on the Kahler potential /C^^-* of the supermul- 
tiplets li and note that for our purposes we can neglect such terms, with the exception 
of the right-handed sector, where terms up to the second order in t have to be taken into 
account. 

At the first order in u we have: 

ka - kA - k' - k' - 

;C« = f^{Mll)s) + y^{MII)a) + f{VT{ll)s) + f{^Till)A) + h.c. (25) 

(■ ■ ■) denotes an invariant under A^, while (■ ■ ■)' and (■ ■ ■)" stand for the 1' and 1" singlets. 
Finally two ^4 triplets, a and 6, can be combined into the symmetric triplet {ab)s or the 
anti-symmetric one {ab)A- The SU(2) singlet fields l'^ are not affected by the first-order 
correction /C*^^^. 

At the second order in u we have a richer structure: 

^(2) = 4') + 4') , (26) 

with the labels L and R referring to lepton doublets / and singlets l'^, respectively. For 
lepton doublets we find: 

lCf = j2%Wl) (27) 

1=1 / 

where X is the list of Zs-invariant operators, bilinear in the fiavon superfields (ps^r and ^ 
and their conjugates, 

X = {^^, (Pr^ (^T)^ ^T<^r, Vs^s, Ivs, ^s^} > (28) 

and each quantity ki represents a list of parameters since there can be different non- 
equivalent ways of combining Xj with II to form an y44-invariant. There are also obvious 
relations among the coefficients ki to guarantee that IC\^ is real. Note that the sum in eq. 
(27) runs over all bilinears which can couple to form A4-invariants. Whether they lead to 
a non-trivial contribution to the relevant terms depends on the VEVs of the flavons. 

For lepton singlets, we can distinguish a diagonal contribution and a non-diagonal one: 

4'^ = + [4'^]w (29) 

1 ^ 

[4^]'^ = 12^ [(''eUX^)e^^^ + (fc^).(X,)/i>^ + ik^MX^)T^r^] (30) 

/ i=l 
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= E 

i=2 
5 



i=2 

5 



e^r"^ + h.c. 



(31) 



i=2 



AJ 



-(Xi)'6'FAr/i'=r'^ + /i.e. 



Due to the structure of the flavon VEVs only the term with 2 = 5 gives a non-vanishing 
contribution in the sum in eq. (31). In order to generate a set of super symmetry breaking 
soft mass terms, the quantities: 



ki = {ko, {k^)i, CK)i, ku,d, kpN, ks, kA, k's, k'^, kj, {kl)u {k'pi, {kl)u {kl^)u {kl,)h (^lir)^ 



are treated as superfields with a constant 9'^ 9 component: 



(32) 



(33) 



where pi and qi do not depend on the Grassmann variables |^ The quantities pi and qi are 
parameters with absolute values of order one. In particular, it is not restrictive to choose 
(see eq. (pil)): 



ko 



1 + go e''9 

^'sUSY 
^2 



{ko)i 



1 + {ql), 9^9 ml^sY 



(34) 

ku,d = 1 + qu,d S'^S^mljjgY , kpN = 1 + qPN (^'^d^rnl^gy . (35) 

When the fiavon fields acquire a VEV according to the pattern shown in eqs. (g and 
eq. (12), the Kahler potential /C gives rise to non-canonical kinetic terms for lepton and 
slepton fields of the following formj^ 



(36) 



The matrices K and K'^ are given by: 



1 + 2ti u 



t4 U 



K 



tiV? 1 — (ti + t2) U 



(ti - ta) u 



(37) 



®In principle we could allow for a more general expansion in the Grassmann variables 9 and/or 9, 
including also terms proportional to 9^ and to 9 . These terms can be absorbed in our parametrization 
after a suitable redefinition of the parameters. 

^We adopt two-component spinor notation, so for example e (e'^) denotes the left-handed (right-handed) 
component of the electron field. For instance, in terms of the four-component spinor ip'^ = (e e'^), the 
bilinears ea'^e and e'^cr'^e'^ correspond to ip^j'^PL^'e and ^^l^^R^e [Pl.r. — ^(lT7^)] respectively. We take 



cr'^ = (l,a), CT^ = (1,-ct), a^'" = \{a^'a''-a''a^'), a 



fll' — 1 



{at^a'^-a^at') and = diag(+l, -1, -1, -1), 



where a = {cr^,a'^,a^) are the 2x2 Pauli matrices. Here the four-component matrix 7^ is in the chiral 
basis, where the 2x2 blocks along the diagonal vanish, the upper-right block is given by cr^ and the 
lower-left block is equal to a^. 
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l + t'iu^ + t'l tl uH tl uH^ 

K''= \ ¥i uH l + t^^u'^ + t'l e tl uH I . (3^ 

F5 uH'^ ¥e uH l + tlu^ + 1% t' 

The coefficients ti^2, ^123 t'123 are real, while the remaining coefficients are complex. 
As one can see, the corrections to the kinetic terms which render these non-canonical 
are small, at most at order 0{u) and O(t^). They have to be taken into account when 
calculating lepton and slepton masses. The coefficients ti and tl are linearly related to the 



parameters pi introduced before, see eq. (33 ). Such a relation is not particularly significant 
and in the rest of this paper we will treat ti and t^ as input parameters, with absolute values 
of order one. 

Notice that we have neglected possible sub-leading contributions to the kinetic terms of 
Hu di and to the fiavons themselves and 9fni which have no impact on the present analysis. 



Note further that through the expansion of ku^ as given in eq. (35) also the soft masses 
'm?H^i for the Higgs doublets Hu and Hd are generated of size qu,d^susY- Similarly, the 
expansion of kpN gives rise to a soft mass term for 9fn of size qpN'm'susY- 

2.2.2 Superpotential 

We continue with the discussion of the superpotential: 

w = wi + Wy + Wd + Wh + ... (39) 
There is a part responsible for the charged lepton masses: 

wi=wf^ +wf^ + ... (40) 
The leading term in the u expansion is: 



= ^(^'pNe^Hd (ipTl) + ^ep^^i^Hd i^Tiy + ^T'^Hd i^Tl)" . (41) 

If ^ 



At the next order in u we find: 

= ^GlNe'Hd (^|/) + ^Op^fi^Hd i^ll)' + ^r^Hd (^|/)" . (42) 



'V 'V 'V 

To generate both the Yukawa interactions and the soft mass contribution ij^eRL)^i 
regard the quantities Xf and x'^ as constant superfields, of the type: 

= yf- ZfO'^msusY , x'f = y'f- z'fO'^TnsusY if = e, /i, r) , (43) 



where the coefficients yf, Zf, y'j and z'j have absolute values of order one. From eqs. (41) 



and (42), after the breaking of the fiavour and the electroweak symmetries, we find the 



following mass matrix for the charged leptons: 

yet'^u + {y'^ + d ye)t^u^ c yet^v^ c yet^v^ \ ^^^g^ 

mi= \ c y^tv? y^tu + {y[^ + c' y^)tu^ c y^^tu^ — -j=- 

c y-rV^ c yrU^ yrU + (y; + d yr)u^ I 

(44) 
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Note that the matrix mi is shown in the basis in which the kinetic terms for the lepton 
fields are given by eqs. ( [36|37[38[ ) . At the LO in the u expansion, rrii matches the matrix 



shown in eq. (13). The off-diagonal elements of mi, all proportional to c, originate from the 
sub- leading contributions to the VEV of the ipT multiplet. Similarly, the superpotential 
giving rise to neutrino masses can be expanded as: 



= wW+w;(2) + ... (45) 



with the LO terms: 



w^J'^ = i^CiHJHJ) + -^i^sHJHJ) , (46) 

AfAL AfAL 

where is the scale at which lepton number is violated. Also in this constant 
superfields: 

Xa,b = ya,b + Za,b0^msUSY ■ (47) 



The terms in w^^ lead to the form of m^, as in eq. (14) with a = UaCa and h = UbCb- The 



next term in the expansion Wu is not relevant for our discussion here, however gives rise to 
deviations of relative order u from TB mixing. The form of w^, responsible for the vacuum 
alignment of the flavon fields, has already been displayed above in eq. (|8]). Finally, the 
term Wh is associated with the parameter: 

Wh = fiHuHd . (48) 

This term explicitly breaks the (continuous) U(l)^ symmetry of the model, while preserving 
i?-parity. The /i term might originate from the Kahler potential of the theory [23], after 
SUSY breaking: 

J J d^9d'^e{XHuHd + h.c.) (49) 

where A is some large scale as, for instance, the Planck scale and X is a chiral superfield, 
whose F component Fx develops a VEV. This gives rise to /i = (Fx) /A. In our model 



we simply assume the existence of the term in eq. (48) in the superpotential and allow 
for an explicit breaking of the U(l)^ symmetry controlled by the parameter /i. The soft 
SUSY breaking term can then arise from the /i term by considering /i as superfield 
+ O^Bfi. A source of the /i term within our model are terms which involve one driving 
field, a certain number of fiavons and the two Higgs doublets Hu and Hd- The lowest order 
term in the superpotential w, allowed by all symmetries of the model, is: 

{^l^T)HuHd/Af. (50) 

This term generates a contribution to the term of the order of msusY x when the driving 
fields acquire VEVs through the inclusion of soft SUSY breaking terms into the fiavon 
potential [16]. The size of such a term is expected to be < 50 GeV for msusY ~ O(lTeV). 
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2.2.3 Slepton masses 

From the Kahler potential and the superpotential we can read off the slepton masses; they 
can be parametrized as follows: 



3 3 



i,j=l i,j=l 

(51) 

where ni?^^^-^^^ and mlj^^ are hermitian matrices and w^eLR — {^IrlV- sneutrino 
sector only the block m^^^ is present. We neglect contributions to the sneutrino masses 
associated to Wy. Each of these blocks receives several contributions: 

tt^Irr = irnlj^R)K + imlj^j^)F + imlj^j^)D + imlj^j^)D,FN ■ (52) 

The contribution to the slepton masses from the SUSY breaking terms in the Kahler 
potential is given by: 

riQ + 2ni u \ 
{'rnlLLlK = {mlLL)K =1 ^4 - (^i + ns) u uqu'^ nilusY 5 

^5 riQ — {rii — 722) u ) 

(53) 



k^lRR)K 



nl v^t 



n^ 



m 



SUSY 



(54) 



n"^ uH 

n^5 uH^ n^e uH 

where the coefficients are complex, except for 720,1,2 and n\ 2 3- p| The coefficients ni and nl 
are linearly related to the parameters qi introduced in eq. (33). Again, such a relation is 



not particularly significant and in the rest of this paper we will treat Ui and ti^ as input 
parameters, with absolute values of order one. 

The SUSY contribution from the F terms is completely negligible for sneutrinos, i.e. 



[m 



uLL)F 



0. For charged sleptons {ml^^)F and {rnlj^j^F read: 



mi 



[m 



eRRJF 



-l\T 



ml 



(55) 



where mi is the charged lepton mass matrix and K, K'^ are the matrices specifying the 



kinetic terms, eqs. (37) and (38). 



The SUSY D term contribution is: 



— - + sin^ 9w 



[m 



uLL)D 



cos 2/5 m\K 
cos 2/3 m^K , 



(56) 



[m 



eRRJD 



-l)sin2 Ow cos 2/3 ml{K' 



c\T 



® Additionally, we assume that no and nl 2 3 are positive in order to have positive definite square-masses, 
to avoid electric-charge breaking minima and further sources of electroweak symmetry breaking. 
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where mz is the Z mass and 6w is the Weinberg angle. Notice again the presence of the 
matrices K and K'^. 

For the right-handed charged leptons we find an additional D term contribution stem- 
ming from the fact that ^fat, and /i"^ are charged under the Froggatt-Nielsen symmetry 
U(l)^^, which we assume to be gauged. The relevant contribution of the U(l)^^ group 
to the scalar potential (through a D term) is: 

1 / „ . r)lC. \^ 

Vd,fn 



1 f ..2 

- [Mpj + gFNQFN-g^Z, 



+ QFNfn'sUSYl^FNl'^ 



(57) 



where stands for the FN charge of the scalar field Zi, and in the second line we have 
displayed only the leading contribution to the terms quadratic in the matter fields. One 
can check that this contribution is of a similar form and size as the one originating from 
the Kahler potential, {fn^RiilK- Thus, we can simply absorb it by redefining the -anyway- 



unknown coefficients n^, i 



1, ...,6 parametrizing {mlj^j^K in eq. (54). 



Concerning the block mlj^j^, this receives three contributions: 



RL 



(m 



eRL)l RL)3 



{51 



The first one originates from the superpotential, eqs. (41) and (42), and is proportional to 
the parameters Zf and z'j: of the decomposition in eq. (43): 



("^eRL)l 



vcosP 
Ai—j=-msusY 



(59) 



with 



Zet^U + (Zg + c' Ze)t^V? C Zet^U^ 

C z^tv? z^tu + {z'^ + c' z^)tu^ 



C Zet^U^ 

C z^tv? 



(60) 



C ZrU^ C Zt-V? ZrU + {z'^ + c' Zr)u^ 

The second one is related to the fact that the auxiliary fields of the fiavon supermultiplet 



ifT acquire non- vanishing VEVs of the form as shown in eq. (21 ), when soft SUSY breaking 



terms are included into the fiavon potential. Formally the contribution can be written as 



dw 



> <9v5t / \ difTdeldlj 
so that the soft mass matrix {mlj^i)2 reads: 

V cos P 

71" 



e^Cj + h.c. 



[m. 



eRL)2 



-^SUSY 



(61) 



(62) 



with 



A2 = C 



cf Vf^tu^ 



Cf yet^u^ 



y^tu + {2y' + dp y^)tu'^ 



Cf VrU^ 



Cf yet^y^ 

Cf y^ltu^ 
yrU + (2?/; + dp 



(63) 
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The third contribution is proportional to the charged lepton mass matrix and is enhanced 
in the large tan fi regime: 



We remark that the origins of the second and the third contribution are quite similar, since 
both arise from the auxiliary component of a superfield: the second one is attributed to 
the auxiliary component of the fiavon superfield (/9t, while the third one originates from 
the auxiliary component associated to the Higgs doublet H^- 



3 The physical basis and its stabihty under renormal- 
ization group running 



In this section we first discuss the results for the slepton masses in the physical basis and 
comment on results found in the literature. In the second part, we give an estimate of the 
renormalization group effects on the slepton masses and show in the leading logarithmic 
approximation that these effects can be neglected or absorbed into the parametrization of 
the soft mass terms. 

3.1 Slepton masses in the physical basis 

All matrices above are given in the basis in which the kinetic terms of the slepton and lepton 
fields arc non-canonical. To derive the physical masses and the unitary transformations 
that enter our computation, we have to go into a basis in which kinetic terms are canonical, 
for both, slepton and lepton, fields. Subsequently, we diagonalize the mass matrix of the 
charged leptons via a biunitary transformation. To avoid flavour- violating gaugino-lepton- 
slepton vertices in this intermediate step, we perform the same transformation on both 
fermion and scalar components of the involved chiral superflelds. This procedure, described 
in detail in appendix B, gives us the physical slepton mass matrices mff.^^n, itt-Irr and 
Th?eRL- The results shown here are obtained under the assumption that all parameters of 
the model are real. The analytical expressions for the slepton mass matrices in the physical 
basis contain the first non- vanishing order in each of the matrix elements. We start with the 
left-left (LL) block. The contribution from the soft breaking terms is common to charged 
sleptons and sneutrinos and reads: 




(64) 




fii — rii — tiUo for i = 1,2, 4, 5,6 



(65) 



The supersymmetric F and D term contributions are given by: 





(66) 
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and 



-- + sm 



Ow COS 2/3 m| X 1 



COS 2/3m|xl , (67) 



with rhi being the mass matrix for the charged leptons in the same basis, i.e. diagonal 
and with canonically normahzed kinetic terms. The supersymmetric D term contributions 
are proportional to the unit matrix. Notice that in the physical basis all SUSY contri- 
butions are diagonal in flavour space. Both, the F and the D term, contributions are 
small compared to that coming from the Kahler potential. The relative suppression is of 
order mJrhi/mgijgY and rn^z/mgugYi respectively, which do not exceed the per cent level 
for typical values of msusY around 1 TeV. Note also that the SUSY part is the only one 
that distinguishes between charged sleptons and sneutrinos. The dominant mass matrix, 
i^eLL)K = {i^IllIki has a structure which is very similar to that of the corresponding 
matrix in the original basis, i.e. all matrix elements in the two bases are of the same order 
in the {u, t) expansion. They only differ for coefficients of order one. 
For rrigjj^ we flnd that {rhlj^j^x is given by: 



K 



m, 



2 c {nl — 713) — u 2 c {nl — 77.3 



2 c {nl — n'. 



rrif. 



-u 



rin 



Trip 



2 c {nl — n'^) — u 2 c {7X2 — ng) 
The supersymmetric terms are: 



m,. 



-u 



rrte 



-u 



\ 



-u 



m 



SUSY 



(68) 



{'>^Irr)f = rhimf and {rhlj^ji)D = 



sin^ 9w cos 2(5 m| x 1 



(69) 



Also in this case the SUSY contributions are diagonal and numerically negligible in most 
of our parameter space. The dominant contribution is thus {Tn1iiji}K- We note that 

)k does not depend on the parameters describing the 



{e,u)LL 



{'^'eRR}K at variance with 
non-canonical kinetic terms. Comparing the size of the entries of {mljiji)K with those of 
("^eiii?)^ we see that the diagonal elements are still of the same order in t and u, whereas all 
off-diagonal elements are enhanced by a factor 1/u. This can be understood in terms of the 
rotation done on the right-handed leptons to diagonalize the charged Icpton mass matrix. 
Such a rotation is characterized by small angles proportional to c ye/y^ tu, c ye/yr t^u and 
c yfji/yr tu in the sectors 12, 13 and 23, respectively. By making the same rotation on the 
corresponding right-handed sleptons, we obtain the off-diagonal terms of rhl^j^. 

Finally, coming to the RL block of the mass matrix for charged sleptons, we find: 



(^eiiL)l = 



me 



yf. 



m-rU 



r, {ZeVn-Zt^ye) „ (Zeyr-Zrye) \ 

2c — nieU 2c nieU 



-m, 



Vt 



rrin-u 



Veyr 



2c 



{Zf,yr - Zry^) 



Zt 



Vt 



msuSY , 



(70) 
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/ nie (cp — c) nieU [cp — c) \ 



msusY , (71) 



{CF — c)m^u {Cp — cjm^u 

\ {cp — c) rrirU {cp — c) rrirU rrir ) 

and 

{mljij^)s = -IJ, tan P mi . (72) 

The matrix ml^j^, which is the sum of these three contributions, does not depend on the 
parameters describing the non-canonical kinetic terms through K and K'^. An important 
feature of (^g^j^)i is that the elements below the diagonal are suppressed by a factor 
u compared to the corresponding elements of (^g/j/,)!, i.e. before the transformations for 
canonical normalization of the kinetic terms and diagonalization of the charged lepton mass 
matrix are applied. However, this does not happen for the second contribution associated 
to the non-vanishing VEVs of the auxiliary components of the supermultiplet (fp so that 
the elements of (?7ig^^)2 are still of the same order in the expansion parameters t and u as 
those of the matrix (^e^i)2. Nevertheless there are cases in which this contribution can be 
suppressed. In the first case the VEVs of the auxiliary fields contained in the supermultiplet 
i^T vanish, i.e. the parameter C, is zero, due to the fact that the soft SUSY breaking terms 
in the flavon potential are (assumed to be) universal, that is equal to the terms of the 
superpotential Wd up to an overall proportionality constant [16]. The second possibility 
arises, if the VEVs of the auxiliary fields can be completely aligned with those of the flavon 
(fp at LO as well as NLO, such that cp becomes equal to c. In both cases the off-diagonal 
elements of (^7ig^^)2 are further suppressed than shown in eq. (71). We emphasize this fact 
here, since it turns out that the suppression of the off-diagonal elements below the diagonal 
as it occurs in the case of is relevant for the actual size of the leading behaviour of 

the normalized branching ratios Rij with respect to the expansion in u. As we shall see in 
section 4, in a general case Rij oc u"^ holds, whereas, if the contribution in eq. ( 71 ) vanishes 
or is also suppressed, Rij is proportional to u"^. The contribution {mlpi^p)^ is diagonal in 
flavour space. Concerning the possible size of this contribution, note that \fi\ tan (3 /msusY 
is the relative magnitude of the non-vanishing elements of {rhlj^^)^ with respect to the 
corresponding ones in (?Tig^j^)i 2. Notice finally that the (31) and (32) element of mlup 
coincide. 

We can compare our results with those found in [24] , where the slepton mass matrices for 
a model possessing the same flavour symmetry were also estimated in a similar framework. 
The main difference between the two setups is that in our model SUSY is a softly broken 
global symmetry, whereas in [24] the model has been embedded into SUGRA. We agree on 
the structure of the matrix (?7igj^j^)x {m\ in the notation of [24]). Concerning the matrix 
{rhl^^K {^"r in [24]), we see that the off-diagonal matrix elements of m|, are all of order 
u^, whereas we find that those of {mljip)K are of order u. Such a discrepancy has a minor 
impact on the estimate of the rates for the radiative transitions, since, as we shall see in 
the next section, the RR block gives a subdominant contribution. It is interesting to note 
that also in the SUGRA context analyzed in [24] the VEVs of the auxiliary components of 
the flavon supermultiplets give rise to {jh\pi)2i ~ rn^nisusY u, which corresponds to our 
("^ei?L)2,2i ~ '^ii'n^susY u, with similar implication on the rate of /i — 67. 
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3.2 Estimate of renormalization group effects 

Here we briefly comment on possible effects of the running from the scale Af ^ at which 
the sfermion masses originate in our effective theory, down to the low energy scale, at which 
the amplitudes of LFV transitions are evaluated. The renormalization group equations for 
the soft mass terms (m^e^^)^^)^, (^ej?H)-'^ and = V2[(m2^^)i + (m^^i)2]/(f cos/?), 
denoting t' = log^Ai /nisusY), are [25]: 

+ 2 (Y}{rnl^^)K% + mlY}% + AlA. 



+ 4(y;(m^iJi.r; + mi/eF; + i,itj^. , (73) 

le^r^l; (^e)^^. = {-lgl-?>gl + mnY,) + MY}Ye))[A:)^^ 

+ 2 0^?iMi + ?>glM2 + 3Tr(F;i,) + 1t{Y}AS^ 

leTT^Ay;., = (^-^^?-3^?2' + 3Tr(f,fJ) + Tr(fen^)^ y;,, 

+ 3 (FeF^Fe' 



where g\^2 are the gauge couplings|^of SU(2)i x U(l)y, Mi 2 the corresponding gaugino 
mass terms, 1^,^ = ^^"^14/ cos /5) are the Yukawa matrices for charged leptons and down 
quarks. Ad = V^[{rhlj^j^) i + (m^^^)2]/(t; cos/3) and: 

S = Tr(mJ^^ + m^^^ - 2m2^^ - (m^^^)^ + (m^^^);^) - m^^ + m^^. 

The matrix {rh'li^i^)K coincides with {rh1j^^)K and has the same evolution. For squarks we 
have introduced soft mass terms analogous to those previously discussed for sleptons. To 
estimate the corrections to the slepton masses induced by the renormalization group evo- 
lution we adopt the leading logarithmic approximation and substitute each of the running 
quantities with their initial conditions at the scale Aj in eqs. (73). In particular, for 



the matrices (?^Li/)ll)^' ('^eRij)^^ ^"^^ Ae, these initial conditions are given in eqs. (65), 
(68) and (70), (71), respectively. In this approximation one easily sees that the largest 
corrections to the matrices {.i^\eu)LL)K a-^id ifhlp^f^K come from electroweak gauge inter- 
actions and are proportional to the identity matrix in flavour space. Due to the negative 

^In the GUT normalization, such that 52 = 9 and gi = y^5/3g' . 
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sign of the dominant contribution these diagonal elements increase by evolving the mass 
matrices from the cutoff scale down to the electroweak scale. This effect is taken into 



account in our numerical study presented in section 5. Each off-diagonal element of 
ifn\ev)LL)K receives at most a relative correction of order: 



u log , (74) 



167r2 \msusY J 

while those to {mlj^^K are even more suppressed, i.e. 



^ log {^^) ■ (75) 
IGtt^ \msusY J 

All such contributions can be safely neglected. 

The matrix A^, gets a first correction by an overall multiplicative factor that can be 
absorbed, for instance, by a common rescaling of the parameters, and a second correction 
of the type ^ Ae + K Y^^ which can be absorbed by the redefinition 
where K and k are constants. In our numerical study these effects are treated in this way. 
Finally, additional corrections to the off-diagonal elements of A^ are negligible. We can 
conclude that the corrections induced by the RG running are either negligible or could be 
absorbed in our parametrization. Thus, in our model the soft mass terms of sleptons are 
completely controlled by the flavour symmetry and by its spontaneous breaking. We recall 
that our model does not contain right-handed neutrinos and ignores the dynamics above 
the scale Ki k, Kj. Notice that the Yukawa couplings of the charged leptons, Y^., remain 
diagonal during the evolution. 

In a see-saw version, one should also include the effects of the running from the cutoff 
scale down to the right-handed neutrino mass scale (s) and the corresponding threshold 
effects. In this case the previous conclusions might change, since we expect for generic 
order one contributions from RG running that they enter the amplitudes of the branching 
ratios at the level of l/(167r^), due to the loop suppression, which is (roughly) equal to a 
contribution of order in the amplitude in our context. 



4 Results in the mass insertion approximation 

We can now evaluate the normalized branching ratios Rij for the LFV transitions — > 67, 
T ^ fi'j and r — *• 67. In this section we estabhsh the leading dependence of the quantities 
Rij on the symmetry breaking parameter u. We recall that in the class of models based on 
the flavour symmetry group A4X Z^xU (1) fn an estimate based on an effective Lagrangian 
approach suggests that Rij generically scales as u^, which can be reconciled with the present 
bound on i?^e only if the scale M of new physics is sufficiently large, above 10 TeV. The 
effective Lagrangian approach also indicates that, under certain conditions, in the SUSY 
case a cancellation might take place, and Rij might scale as a combination of two terms, one 
proportional to and one proportional to {w?- / m'f)u^ , as shown in eq. ([3|, thus allowing 
for a substantially smaller scale of new physics, in the range of (1 ~ 10) TeV. When i?^e is 

^"This effect can also be viewed as a redefinition of the initial parameters riQ and nl 2 3- 
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dominated by a one-loop amplitude with virtual particles of mass msusY, M and msusv 
are roughly related by M = {4:7!-/g)msusY and a given lower bound on M corresponds to 
a lower bound on msusv one order of magnitude smaller. Here we determine the actual 
leading order behaviour of Rij by expressing the result as a power series in the parameter 
u. 



4.1 Analytic results 

In order to do so it is useful to first analyse the predictions in the so-called mass insertion 
(MI) approximation, where we have a full control of the results in its analytic form. A 
more complete discussion based on one-loop results can be found in section 5. For the 
case at hand, the MI approximation consists in expanding the amplitudes in powers of the 
off-diagonal elements of the slepton mass matrices, normalized to their average mass. From 
the expression of the mass matrices of the previous section we see that in our case such 
an expansion amounts to an expansion in the parameters u and t, which we can directly 
compare with eq. (|3]) . A common value in the diagonal entries of both LL and RR blocks 
is assumed and we consequently set uq = n1 = n2 = n'^ = 1 and also fii = f;,2 = in 
this section, so that the average mass becomes msusv- On the contrary, no assumptions 



have been made for the trilinear soft terms, which keep the expression as in eqs. (70 



72). Concerning chargino and neutralino mass matrices, they carry a dependence on the 
vector boson masses mw,z through off-diagonal matrix elements. Such a dependence is not 
neglected in this approximation, but only the leading order term of an expansion in mw,z 
over the relevant SUSY mass combination is kept. At the same time, to be consistent, we 
have to neglect the supersymmetric contributions of m^i^i^ and ^eLL ^^"^ therefore m^i^j^ 
and ml^^ coincide. Using these simplifications, the ratios Rij can be expressed as: 

^.• = 7^^(i^!p+Kn • (76) 

'^piii'SUSV 

At the LO, the amplitudes A^l and are given by: 

_ /X \ , f^SUSYfr- N 

~ O.LL\0ij)LL + O-RL [OijJRL 

m^ 

= aRR{Sij)RR + aiR ^^^^^ {Sij)LR (77) 

mi 

where acc' {C, C = L, R) are dimensionless functions of the ratios Mi^2/fnsusv, l^/i^susv 
and of tan^vy- Their typical size is one tenth of / {IQn'^) , g being the SU(2)i gauge 
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coupling constant. In our conventions their explicit expression is given by: 



(III = 



167r2 



hn{a2) + flc{a2) + ^2 _ ^2 



f2n{a2,b) + /2c(a2,&)) 



2„ / ^ / X Ml// tan/? 
+ tan Ow /in(ai) 7-3 ^/2n(ai, 0) - Mi 



(iRR — 
dLR = 



Mi^ - //2 

tan e'w 2/2„(ai) 
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where ai,2 = Ml^/ml^sY^ b = fJ-^/mljjgY and fi(c,n)ix,y) = fi(c,n)ix) - fi{c,n){y)- The 
functions fin{x) and fic{x), slightly different from those in [26], are given by: 



fin{x 

f2n{x 

flc{x 
f2c{x 



-17x^ + 9x'^ + 9x-l + 6x^{x + 3) log x) /(24(1 - xf) 
-5x^ + 4x + 1 + 2x{x + 2) log x) /(4(1 - x)^) 
l + 9x-9x^-x^ + 6x{x + 1) logx)/(2(l - xf) 
-x^ -9x^ + 9x + l + 6x{x + 1) logx)/(6(l - xf) 
-x^-Ax + 5 + 2(2x + 1) log x)/(2(l - x)^) . 
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Notice that ace do neither depend on u nor on the fermion masses rriij. Finally, {Sij)cc' 
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Table 2: Coefficients w'?^ characterizing the transition amphtudes for /j, ej, t and t iij, in 

the MI approximation in which no and are set to one and ni,2 to zero so that vanish. 

parametrize the Mis and are defined as: 



m 
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Prom the mass matrices of the previous section, we find (j < i): 

n^susY ^ ' 
m msusY 



(81) 
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where for the mass insertion {5ij)RL we have also displayed the NLO contributions, in order 
to better compare our results with those of the effective Lagrangian approach. The explicit 



expression for the LO coefficients wf^^' are listed in table 2. Also the NLO coefficients 
yj'RL g^j^g dimensionless combinations of order one parameters. By substituting the mass 



insertions of eq. (81) into the amplitudes r of eq. (77) and by using eq. (76), we get: 
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The behaviour displayed in eq. (82) differs from the one expected on the basis of the 
effective Lagrangian approach in the SUSY case, eq. This is due to the presence 

Assuming w^,^ = we recover what is expected from the 



(0) 



of the term • 



oc w. 



RL 

ij 



effective Lagrangian approach in the SUSY case, whereas when wj^^ does not vanish, the 
LO behaviour matches the prediction of the effective Lagrangian approach in the generic, 
non-supersymmetric case, eq. ([2|. As shown in table 2, the coefficient wfj'" is universal, 
namely it is independent from the flavour indices and it vanishes in two cases: 

i) cp = c, which reflects the alignment of the VEVs of the scalar and auxiliary compo- 
nents of the fiavon supermultiplet ^px, see eqs. (g and (|2l|). 



ii) C = which can be realized by special choices of the soft SUSY breaking terms in 
the fiavon sector, i.e. the assumption of universal soft SUSY breaking terms in the 
fiavon potential. 

In our model none of these possibilities is natural, see [16], and both require a tuning of 
the underlying parameters. If wfj^ = 0, the result expected from the effective Lagrangian 
approach in the SUSY case is obtained in a non-trivial way. Indeed, it is a consequence 
of a cancellation taking place when going from the Lagrangian to the physical basis. In 
particular, for w^^ = the MI {5ij)RL scales as rriiU^ /uisusy and not as rriiu/msusY as 



we might naively guess by looking at eq. (60). As a consequence RfF^^ scales as u'^ and 



not as u"^ for rrij = 0. 

In the general case when wjj^ is non-vanishing, the dominant contribution to RfJ^^"^ 
regarding the expansion in u is flavour independent and, at the LO in the u expansion, we 



predict R^ 



pe 



R. 



variance with the predictions of most of the other models, 
where, for instance, R^e/Rrii. can be much smaller than one [1,27,28]. If w^^ is non- 
vanishing, it is interesting to analyze the relative weight of the leading and sub-leading 
contributions to Rij. For this purpose we list in table |3] the expressions and the numerical 
values of the functions acci in the limit yU = Mi^2 = ^susy- As one can see, in this limit 
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Table 3: Coefficients ace characterizing the transition amplitudes for fj, e^, t and r —>■ fij, 

in the Mf approximation in which uq and are set to one and hi,2 to zero and by talcing the limit 
the /i — All 2 — ^susY- Numerical values are given in units of g^/{192TT^) and using sin^ 6y\,' — 0.23, 
{zi/ui + C) = 1; Eind tan f3 ~ 2 ^ 15. 



the dominant coefficient is ull, which is larger than a^j^ = ulr by a factor 7^54, and larger 
than aRR by a factor — (4 -h 13), depending on tan/3 = 2 15. Assuming coefficients wfj'^' 
of order one in eqs. (83), we see that the most important contributions in the amplitudes 
for the considered processes are aji^u and anv?. The ratio between the sub- leading and 
the leading one is {aLi/ aRL)u ~ (7 ^ 54)^. When u is close to its lower bound, which in 
our model requires a small value of tan /3, the leading contribution clearly dominates over 
the sub-leading one. However, for u close to 0.05, which allows to consider larger values of 
tan/3 K, 15, the non-leading contribution can be as large as the leading one and can even 
dominate over it. The transition between the two regimes occurs towards larger values of 
u. 

The numerical dominance of the coefficient an has also another consequence: for 
vanishing w^^ , Rij is dominated by the contributions of aLLwfj'", whose values are not 
universal, but expected to be of the same order of magnitude for all channels. Thus even 
when wfj^ = 0, we predict R^e ~ Rrfi ~ Rre- Even if to illustrate our results explicitly 
we have taken the special limit /i = Mi 2 = fnsusY, a numerical study confirms that they 
remain approximately correct when more generic regions of the parameter space of the 
model are considered. 



4.2 Failure of the effective Lagrangian approach 

Why the results of the effective Lagrangian approach in the SUSY case do not apply in the 
model under consideration? It is worth to summarize the main assumptions underlying 
the effective Lagrangian approach in the SUSY case in [4]: 

1. The only sources of chirality ffip are either fermion masses or sfermion masses of 
RL-type. 
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2. Both, fermion masses and sfermion masses of RL-type, up to the order u^, are dom- 
inated by the insertions of cpT or (p"^ in the relevant operators. 



As a consequence also the operators describing dipole transitions, which have the same 
chiral structure as the charged lepton mass terms, are dominated by the insertions of ifT 
or (p"^, up to the order u"^. For instance one such operator is: 

' ^^^e^/7.a,.F-f^A+- (84) 



^SUSY . 

and similarly for /x"^ and t'^. Dots denote additional contributions such as those coming 
from the insertion of y?^. Here e is the electric charge. 

Indeed, in our model condition 2 is violated. The explicit SUSY model considered here 
contains another set of fields which was not present in the effective Lagrangian framework, 
namely the driving fields ipQ. The driving fields represent an important ingredient of our 
model, since they are directly responsible for the vacuum alignment of the flavon fields. 
In the effective Lagrangian approach the alignment was postulated, without referring to a 
specific dynamical mechanism to generate it, and the driving fields were not included among 
the relevant degrees of freedom. In our model (p^ has no direct coupling to matter fields. 
Such a coupling arises indirectly through the mediation of the flavon ifT, which is coupled 
to matter via interactions suppressed by 1/Aj and to via interactions proportional to 
e.g. a large scale M^: 

w = xfi^e'^H, [^^l^ + Mt{^I^t) + ... (85) 

Mt is a mass scale of the order of the VEV of (pT-, Mt ^ u Af and dots denote further 
contributions. It can be shown that the VEV of <^o is proportional to msusY [16]. When 
the auxiliary fields of (pT are eliminated slepton masses of RL-type receive a contribution 
from the insertion of • 

"FN ~CTT T^j yor 



ye-j^~e^H,^MT^lj+... (86) 

Similar terms for fjif^ and are also present. Therefore dipole operators with the same 
type of insertions are expected in the low-energy limit. Up to loop factors: 

-^'i^e'H.a^.F'"- (Mrfl) + ... (87) 

''''SUSY ^^f \ ^^f / 

Since has a dominant VEV of order msusY and Mt/Aj ^ u, this operator is similar in 



size to the operator of eq. (84). A similar contribution arises, if the coupling of two flavons 
ipT to the driving field (pQ is taken into account instead of the term Mt{'Pq'Pt)- Operators 
of this type were not included in the effective Lagrangian approach, simply because the 



driving field ip^ was absent. The set of operators of type as in eq. (87) and eq. (84) 
are aligned in flavour space only to the LO. At the NLO such an alignment fails and this 
produces a no n- vanishing term wfj^. Notice that the insertion of (Pq both in the slepton 
mass and in the dipole operators occurs through a combination of the type MrpTp^/Af. If 
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we did not account for the large coupling of order M^, the mere insertion of Tp'^ /Kf in 
the mass and in the dipole operators would largely underestimate the effect. It would be 
interesting to investigate the low-energy structure of dipole operators in a complete SUSY 
model in which the driving fields are absent and the alignment of flavon vacua is realized 
via some alternative mechanism. This would allow to verify whether the failure of the 
effective Lagrangian approach is entirely due to the presence of the driving fields or other 
subtleties occur when performing the low-energy limit. 



5 Numerical analysis 

In this section we perform a numerical study of the normalized branching ratios Rij and 
of the deviation 5a^ of the anomalous magnetic moment of the muon from the SM value. 
We use the full one-loop results for the branching ratios of the radiative decays as well as 
for 5a^. These can be found in [29-32] and are displayed in appendix C for convenience. 



5.1 Framework 

As discussed in the preceding sections, in our model the flavour symmetry x Z3 x 
U(l)^^ X U(l)j:j constrains not only the mass matrices of leptons, but also those of 
sfermions. These are given at the high energy scale A/ A^,, which we assume to be 
close to 10^^ GeV, the SUSY grand unification scale. The flavour symmetry does not fix 
the soft SUSY mass scale rususY- It also does not constrain the parameters involved in 
the gaugino as well as the Higgs(ino) sector. These are fixed by our choice of a SUGRA 
framework in which mgusY is the common soft mass scale for all scalar particles (in the 
literature usually denoted as mo) and mi/2 the common mass scale of the gauginos. Thus, 
at the scale A j ^ Aj;^ we have 

Mi(Ai) = M2(Ai) = mi/2 . (88) 

Effects of RG running lead at low energies (at the scale mw of the W mass) to the following 
masses for gauginos 

M,{m^) ^ ^^MMl) MMw) - ^^^M2(A,) , (89) 

where ctj = gf/4:7i {i = 1,2) and according to gauge coupling unification at Aj ^ A^, 
Q;i(Ai) = Q;2(Ai) ~ 1/25. Concerning the effects of the RG running on the soft mass terms, 
as we have seen in section 3.2 these are small or can be absorbed into our parametrization 
of the soft mass terms. Thus, in the contributions (^7ig^^)i_2 to the RL block we take 
rrisusY as input parameter. Nevertheless, we explicitly take into account the RG effect on 
the average mass scale of the LL block, m|, and in the RR block, mj^, 



mlimw) ^ rnl{AL) + 0.5M|(Ai) + 0.04M2(Ai) - ml^gy + OMml/^ , 
m|(m^) ~ m|(Ai) + 0.15Mf{AL) ~ ml^gy + 0.15m2^2 • 
^"'^The gluino mass parameter M3 is not relevant in our analysis. 



(90) 
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The parameter /i is fixed through the requirement of correct electroweak symmetry breaking 

, ,9 ^%.~^'h tan^/3 1 „ ^ ^ 

I"! = Un^;il 2'"^ ■ 

Since in the SUGRA framework the soft Higgs mass parameters are also given by msusY 



at the high energy scale, mj^^(Aj^) = mH^{J^L) = ''^susy ' relation in eq. (91) reads at 



low energies 



m| _2 l + 0.5tan2/5 2 0.5 + 3.5 tan^ /? 
tan2/?-l +""1/2 tan2/5 



so that /i is determined by msusY, ?^i/2 and tan/3 up to its sign. We recall that in our 
model the low energy parameter tan P is related to the size of the expansion parameter u, 



the mass of the r lepton and the r Yukawa coupling y^-, as shown in eq. (18). For this 
reason, requiring 1/3 < \yr\ < 3 constrains tan/3 to lie in the range 2 < tan/? < 15. As 
already commented, the lower bound tan /3 = 2 is almost excluded experimentally, since 
such low values of tan (3 usually lead to a mass for the lightest Higgs below the LEP2 bound 
of 114.4 GeV [33]. Q 

In our numerical analysis the parameters are the following: the two independent mass 
scales msusY and mi/2, the sign of the parameter /i and the parameters of the slepton 
mass matrices shown in section 3.1 in the physical basis. We recall that the results of 
section 3.1 have been obtained under the assumption that the parameters are real and 
we keep working under the same assumption here. We also assume that the parameters 
on the diagonal of the slepton mass matrices ifn\ev)LL)K and ('Tig^^)x, ""-o and n^ 2,3' 
positive in order to favour positive definite square-masses, to avoid electric-charge breaking 
minima and further sources of electroweak symmetry breaking. The absolute value of the 
0{1) parameters is varied between 1/2 and 2. We will chose some representative values for 
u in the allowed range 0.007 ^ m ^ 0.05. The other expansion parameter t is fixed to be 
0.05. In the analysis of the normalized branching ratios Rij we fix tan/3 and u and then 
we derive the Yukawa couplings ?/e, Vfj. and Ur- When discussing the anomalous magnetic 
moment of the muon instead we vary Ut between 1/3 and 3 and calculate tan/3 by using 



eq. (18). Having determined tan/3, the Yukawa couplings Ue and can be computed. 

The allowed region of the parameter space is determined by performing several tests. 
We check whether the mass of the lightest chargino is above 100 GeV [18], whether the 
lightest neutralino is lighter than the lightest charged slepton, whether the lower bounds 
for the charged slepton masses are obeyed [18] and whether the masses of all sleptons are 
positive. The constraint on the mass of the lightest chargino implies a lower bound on mi/2 
which slightly depends on the sign of /i. In our plots for Rij we also show the results for 
points of the parameter space that do not respect the chargino mass bound. For low values 
of rrisusY, e.g. msusY = 100 GeV, the requirement that the lightest neutralino is lighter 
than the lightest charged slepton is equivalent to the requirement that the parameters in 
the diagonal entries of the slepton mass matrices u)ll) k and {rhl^j^)K are larger than 
one. For larger values of msusY, e.g. msusY = 1000 GeV, this requirement does not affect 



-"^^This bound assumes that the Higgs is SM-like. For the case of generic MSSM Higgs the bound is much 
lower, 91.0 GeV [34]. 
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our analysis anymore. We note that masses of charginos and neutralinos are essentially 
independent from the 0{1) parameters of the slepton mass matrices and thus their masses 
fulfill with very good accuracy (better for larger mi/2) 

Mjo ^ 0.4mi/2 , Mjo ^ M~- ^ 0.8mi/2 , M^o ^ M^o ^ M~- ^ |^| . (93) 

For the slepton masses we find certain ranges which depend on our choice of the 0{1) 
parameters. 

5.2 Results for radiative leptonic decays 

We first discuss the results for the branching ratio of the decay ^ 67. This branching 
ratio is severely constrained by the result of the MEGA experiment [35] 

R^e ~ BR{fi 67) < 1.2 X 10"^^ (94) 

and will be even more constrained by the on-going MEG experiment [36] which will probe 
the regime 

R^e ~ BRifi ^ 67) > 10-1=^ . (95) 

We explore the parameter space of the model by considering two different values of the 
expansion parameter u, u = 0.01 and u = 0.05, two different values of tan /3, tan (5 = 2 and 
tan/3 = 15, as well as two different values of the mass scale msusY, f^susY = 100 GeV 
and msusY = 1000 GeV. We show our results in scatter plots in Figure [T] choosing mi/2 
to be mi/2 < 1000 GeV. All plots shown in Figure [l] are generated for /i > 0. 

As one can see from Figure [T]^a), for very low tan P = 2, small u = 0.01, small msusY = 
100 GeV the experimental upper limit from the MEGA experiment on BR{fi —>■ 67) can be 
passed in almost all parameter space of our model for values of mi/2 as small as 450 GeV. 
For msusY = 100 GeV and mi/2 = 450 GeV the sparticle masses are rather light: the 
lightest neutralino has a mass of 175 GeV, the lightest chargino of 350 GeV, the masses of 
the right-handed (charged) sleptons vary between 175 and 285 GeV and the masses of the 
left-handed sleptons are in the range (250 500) GeV. Thus, especially the right-handed 
sleptons are expected to be detected at LHC. In a model also including quarks (and hence 
squarks) we find for the squarks that they can have masses > 700 GeV and gluinos with 
masses of about 1000 GeV, all accessible at LHC. To pass the prospective bound coming 
from the MEG experiment in a sizable portion of parameter space of our model mi/2 has to 
be chosen larger, mi/2 > 600 GeV. Then, however, the masses of the sleptons might only 
be detected at LHC in case of right-handed sleptons. As one can see, values of mi/2 < 155 
GeV are excluded due to the constraint on the lightest chargino mass. Studying the same 
value of tan P and u, but taking tususy to be as large as 1000 GeV, we can see from Figure 
|l|b) that now the bound set by the MEGA experiment on BR{^ 67) is respected in 
the whole parameter space of our model for all values of mi/2- Also the foreseen limit of 
the MEG experiment can only exclude a smaller part of the parameter space of the model 
for all values of mi/2- In this setup, the prospects for detecting SUSY particles at LHC 
are the best for gauginos due to the possible low value of mi/2- The slepton masses are 
expected to be roughly msusY and thus too large to allow for a detection at LHC. 
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Increasing the value of the expansion parameter u from u = 0.01 to u = 0.05, as 
done in Figure [T|c) and[l|^d), increases also the branching ratio of the decay /i — 67 by 
approximately two orders of magnitude, since the different contributions to the branching 
ratio scale at least with m^, as analyzed in section 4. For this reason for low values of 
^susY = 100 GeV, mi/2 has to take values mi/2 ^ 600 GeV in order for the result of 
BR{fi 67) to be compatible with the MEGA bound at least in some portion of the 
parameter space of our model. For the point {msusY, ^1/2) = (100 GeV, 600 GeV) the 
sparticle spectrum is characterized as follows: the lightest neutralino has a mass of 240 
GeV, the lightest chargino of 470 GeV, right-handed sleptons between 250 and 350 GeV 
and left-handed sleptons generally above 300 GeV. For this reason there still exists the 
possibility to detect right-handed sleptons at LHC. Concerning gluinos and squarks these 
are expected to have masses between 1000 and 1500 GeV so that they also can be detected 
at LHC. As one can see from Figure ^c) the on-going MEG experiment can probe nearly 
the whole parameter space of the model for tan/3 = 2, u = 0.05 and msusv = 100 GeV 
for values of mi/2 ^ 1000 GeV. Increasing the parameter msusY to 1000 GeV shows that 
applying the existing bound on BR{^ 67) of 1.2 x 10^^^ cannot exclude small values of 
mi/2- The situation changes, if the expected bound from the MEG experiment is employed, 
because then values of mi/2 smaller than 1000 GeV become disfavoured. 

Finally, we show in Figure [T]^e) and [T]^f) the results obtained for tan/3 = 15. We 
remind that this value is the largest possible one of tan f3 in our model. Requiring that the 
r Yukawa coupling does not become too large entails that tan = 15 fixes the expansion 
parameter u to take a value close to its upper limit, u = 0.05. The value of BR{fi — > 67) 
is thus enhanced through tan /3 as well as u. This is clearly shown in Figure [l]^e) and 
[T|f), because for a low value of msusv = 100 GeV already the MEGA bound practically 
excludes almost the whole parameter space of our model for all values of mi/2 < 1000 
GeV. Increasing the mass parameter msusv to 1000 GeV slightly improves the situation, 
because now there exists a marginal probability to pass the MEGA bound. Again, however, 
the MEG experiment can probe all parameter space of our model for mi/2 ^ 1000 GeV. 
Thus, for msusv ^ 1000 GeV and mi/2 < 1000 GeV the parameter space of our model is 
already severely constrained for moderate values of tan/? which entail large u ~ 0.05 by 
the bound coming from the MEGA collaboration, but surely will be conclusively probed 
by the MEG experiment. Choosing /i < hardly affects the results presented here apart 
from slightly decreasing the lower bound on mi/2 coming from the chargino mass bound. 
Thus, all statements made also apply for < 0. 

In summary, the current bound on BR{fi 67) prefers regions in the parameter space 
of our model with small u or small tan (3, as long as the SUGRA mass parameters should 
be chosen smaller than 1000 GeV. The foreseen MEG bound strongly favours regions in 
which u is small for msusv and mi/2 being not too large. The fact that smaller values of u 
are preferred has consequences also for the expectations of the detection prospects for the 
reactor mixing angle 6'i3, because this angle scales with u: it might thus not be possible to 
detect 6*13 with the reactor and neutrino beam experiments under preparation [37,38]. 

Concerning the radiative r decays, t —>■ fi'-f and r ej, the result found in the 
MI approximation that the branching ratios of these decays are of the same order of 
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(e) tan /3 = 15, u = 0.05 and msusv = 100 GeV. (f) tan/3 = 15, w = 0.05 and msuSY = 1000 GeV. 

Figure 1: Scatter plots of BR{ii ej) as a function of mi/2, for different values of tan/3, u and msusv- 
The red (dark gray) points correspond to points in which the mass of the lightest chargino is below the 
limit coming from direct searches. The horizontal lines show the current MEGA bound (continuous line) 
and the prospective MEG bound (dashed hue). 



magnitude as BR{fi — > ej) is essentially confirmed in a numerical analysis. Due to the 
random parameters differences up to two orders of magnitude are expected and found, 
especially for the case of larger tan/5. However, it is still highly improbable that the 
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decays t —>■ jj,^ and r — 67 could be detected at a SuperB factory, assuming a prospective 
limit of BR{t 7x7), BR{t 67) > lO'^ [39]. 

5.3 Results for anomalous magnetic moment of the muon 

As is well known, the value found for the anomalous magnetic moment of the muon [40] 

aJ^-P = 116592080(63) x 10"" (96) 

shows a 3.4 0" deviation 

= aj^^ - a^*^ = +302(88) x 10"^^ (97) 

from the value expected in the SM [41] 

a^*^ = 116591778(61) x 10"^^ . (98) 

Similar discrepancies have been reported in [42]. Thus, it might be interesting to consider 
the case in which this deviation is attributed to the presence of SUSY particles with masses 
of a few hundred GeV. The one-loop contribution to the anomalous magnetic moment of the 
muon in supersymmetric extensions of the SM has been studied by several authors [29]. We 




Figure 2: Scatter plots in the plane BR{fj, ej) — (5a^, for values of m = 0.01, 0.05. The value of tan/3 
is fixed through the relation with the r Yukawa coupling, which lies in the interval [1/3,3]. The values of 
msusY and of mi/2 are chosen between 10 and 300 GeV for the left panel and between 10 and 1000 GeV 
in the right one. The horizontal lines correspond to the MEGA (continuous line) and the MEG bounds 
(dashed line) ; the vertical lines correspond to the measurements on Ja^ : the continuous one is the best fit 
value and the dashed ones correspond to the 3(t boundaries. 



study the compatibility between the requirement that 6a^ is explained by the exchange of 
relatively light SUSY particles and the experimental upper limit on BR{fi 67) coming 
from the MEGA experiment. We choose again two different values of u, u = 0.01 and 
u = 0.05. To better explore the parameter space we vary the r Yukawa coupling between 
1/3 and 3 and fix the value of tan f3 through the relation given in eq. ( 18 ). As a consequence 
in the plot for u = 0.01, see Figure [2|a), 2 < tan/5 < 3 holds and for u = 0.05 tan/3 
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takes values 2 < tan f3 < 15. Similarly, msusY and mi/2 are chosen to lie in intervals 
[10 GeV, 300 GeV] and [10 GeV, 1000 GeV] for u = 0.01 and u = 0.05, respectively. The 
different choice of intervals is due to the fact that values of a few hundred GeV for msusY 
and mi/2 are disfavoured by the existing limit on BR{fi — > ej) when u = 0.05. As one 
can clearly see from Figure [2| in almost the whole parameter space of our model it is not 
natural to reproduce the observed deviation of the muon anomalous magnetic moment and 
at the same time to respect the existing bound on the branching ratio of /i — 67. 

This kind of incompatibility is well known in supersymmetric theories, because the 
explanation of the 3Aa discrepancy necessitates small values of msusY ssid mi/2 and 
larger values of tan/5, which in turn enhance the branching ratio of the radiative LFV 
decays. Thus, we have to conclude that either there exist further sources of contributions 
to the anomalous magnetic moment of the muon beyond those present in our model, or - 
as is also discussed in the literature - the theoretical value af,^'^ found in the SM is closer 
to a^^^ so that the eventual discrepancy becomes less than 100 x 10~^^, a value which 
could well be explained in our model. 

6 Conclusion 

While awaiting the start of LHC and the first exploration of the TeV scale, the solution 
to the hierarchy problem offered by SUSY still represents a very appealing option, with 
several interesting consequences, such as a successful gauge coupling unification, a viable 
particle candidate for dark matter and a possible explanation of the observed discrepancy 
in the anomalous magnetic moment of the muon. Nevertheless the existence of new states 
carrying flavour indices at the TeV scale is a serious problem in model building, given the 
success of the SM in describing all known phenomena involving flavour transitions. The 
lepton sector makes no exception. The flavour conversion observed in neutrino oscillations 
might be related to potentially large effects in rare decays of the charged leptons. A naive 
dimensional estimate of these effects, assuming new physics at the TeV scale, requires a 
large suppression in the relevant amplitudes, to pass the existing bounds. A possible expla- 
nation for such a suppression is the presence of a flavour symmetry, which is independently 
motivated by the hierarchy in the charged lepton mass spectrum and by the nearly TB 
mixing pattern of the lepton mixing matrix. Here we have considered a SUSY model in- 
variant under the flavour symmetry group A^x Z3 x U{1)fn, originally proposed in order 
to describe lepton masses and mixing angles. We have extended the original model by in- 
cluding SUSY breaking terms consistent with all symmetry requirements. Our model is an 
effective theory, valid at energy scales below a cutoff Aj ~ A^, where we have derived the 
spectrum of SUSY particles, in the slepton sector, under the assumption that the SUSY 
breaking scale is larger than Af. It provides an example of a model in which the slepton 
mass matrices at the scale Aj are not universal. Left-handed sleptons are approximately 
universal, with a small departure from universality controlled by u, the flavour symmetry 
breaking parameter. Right-handed sleptons have soft masses of the same order, but the 
relative difference among them is expected to be of order one. Off-diagonal elements in 
both sectors, as well as in the RL block, are small and have a characteristic pattern in 
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terms of powers of u. This structure is maintained by the effects coming from the RG 
running from A/ ^ down to the clcctrowcak scale. The symmetry breaking parameter 
u hes in a restricted range around few per cent and has a size similar to the reactor mixing 
angle ^^13. We have exploited the knowledge of the slepton mass matrices to compute the 
normalized branching ratios Rij for the transitions — > 67, r — > /X7 and r — 67. At vari- 
ance with other models based on flavour symmetries we found R^e ~ it^r/* ~ Rre and, given 
the present limit on R^ei these rare r decays are practically unobservable in our model. 
On a more theoretical side, the scaling R^j oc v?, found in the MI approximation, violates 
an expectation based on an effective Lagrangian approach, which suggested R^j cx in 
the limit of massless final charged lepton. We have identified the source of such a violation 
in a single, fiavour independent, contribution of the RL block of the slepton mass matrix. 
Such a contribution originates from the VEVs of the auxiliary components of the fiavon 
supermultiplets. We have classified the conditions under which this universal contribution 
is absent. 

In a numerical analysis of R^^ we found that already the current bound from the 
MEGA experiment requires the parameter u to be small or tan (3 to be small for SUSY 
mass parameters msusY and mi/2 below 1000 GeV, to guarantee detection of sparticles at 
LHC. Applying the prospective MEG bound tightens the parameter space of our model 
even more to small u and tan/? or requires mass parameters msusY and mi/2 above 1000 
GeV. Furthermore, we showed that the deviation of the experimentally observed value of 
the magnetic moment of the muon from the SM one cannot be naturally explained in our 
framework, for BR[ii — > 67) below the current bound. The maximal value of in our 
model is around 100 x 10"^^ for BR{n 67) < 10"^^ 
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Appendix A The group A4 

The group /I4 is generated by two elements S and T obeying the relations [14]: 

52 ^ (^7^)3 ^T^^l . (99) 

It has three independent one-dimensional representations, 1, 1' and 1" and one three- 
dimensional representation 3. The one-dimensional representations are given by: 

1 S^l T=l , 



The three-dimensional representation, in a basis where the generator T is diagonal, is given 
by: 



r= 0a;2 , S^-{2-12\ . (101) 





The multiplication rule for triplet representations is the following: 

3 X 3 = 1 + 1' + 1" + 35 + 3a (102) 

If we denote by: 

0^(01,02,03) , b^{bi,b2,h) (103) 

two triplets, the singlets contained in their product are given by: 

1 = (ah) = (ai6i + 02^3 + ^3^2) 

1' = (ab)' = (0363 + ai&2 + a2&i) (104) 
1" = (ab)" = (0262 + ai&3 + 0361) 

The two triplets can be separated into a symmetric and an antisymmetric part: 

35 = {ab)s = ^(2oi6i - 0,2^3 - «'3^2, 20,363 - 0162 - 02^1, 20262 - 01^3 - 03^1) 

3a = (o6)a = 2(^2^3 - 03^2, ai&2 - 0261, 0361 - 0163) (105) 

Moreover, if c, c' and c" are singlets transforming as 1, 1' and 1", and a = (01,02,03) is 
a triplet, then the products oc, ac' and oc" are triplets explicitly given by (oic, 02C, 03c), 
(03c', oic', 02c') and (02c", 03c", oic"), respectively. Note that due to the choice of complex 
representation matrices for the real representation 3 the conjugate o of o ~ 3 does not 
transform as 3, but rather (oi, 03, 02) transforms as triplet under A4. The reason for this is 
that T* = U'^TU and 5"* = U^SU — S where U is the matrix which exchanges the second 
and third row and column. 
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Appendix B Canonical normalization of the kinetic 
terms and diagonalization of the charged lepton mass 
matrix mi 

We perform the following transformations on the superfields: the kinetic terms are canon- 
ically normalized, i.e. we find the basis in which the matrices K and K'^, shown in eqs. 
( 37p8 ), are unit matrices, and in a second step the charged lepton mass matrix m/, eq. 



(44), is diagonalized. As explained in the main text, we perform these transformations 
not only on the lepton fields, but also on the sleptons in order to ensure that the gaugino- 
lepton-slepton vertices do not violate flavour at this stage. 

To diagonalize the hermitian matrices K and K'^ we apply the unitary transformations 
W and W: 

W^KW = dia.g and = diag . (106) 

Normalizing K and K'^ requires a rescaling of the fields via the real (diagonal) matrices R 
and R^: 

RW^KWR = 1 and R^iW^)^ K'^W^R" = 1 . (107) 

The superfields / and l'^ {l^ = (e^, /i'^, t'^) as above) are expressed as: 

l = WRl' and r = Vr'=i?T' (108) 

so that the kinetic terms for leptons and sleptons are in their the canonical form, i.e. eq. 



(36) becomes: 



The mass matrices for leptons and sleptons in this basis read: 

rmil = r'R^{W^fmiWRl' =r'm\l\ (109) 

and 
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icm\^ I = [c R^i^W^fm^^WR I' , 

lmlJ=l' RW^ mlLWRl' , (110) 
ic r = f' R%W' fmlj^ W^^R" . 

We diagonalize the resulting mass matrix m'l of the charged leptons by the usual bi-unitary 
transformation: 

t/^m^y = diag(me, m^, rrj-r) (HI) 
and arrive at the mass eigenbasis /" and l'^": 

r = Ur" and I' = Vl" . (112) 



^^Note that the matrix m\j^ stands for both, the mass matrix of the left-handed charged sleptons and 
the one of the sneutrinos. 
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Finally, the slepton mass matrices and are given as: 



Iml^ I = f'iV^RW^ ttiIl WRV] I" 



(113) 



I' 



m 



RR 



I' 



As we assume for the actual calculation of rh\j^, ^"ll '^'rr. ^^at all couplings 
involved are real, the matrices W , W^, U and V turn out to be orthogonal instead of 
unitary. Furthermore, we express the small parameter u in terms of t as m = with x < 1 



according to eq. (15) and eq. (19). We then can do the calculation in just one expansion 
parameter t. ^ In the course of the calculation we pose the following requirements: the 
kinetic terms are canonically normalized up to and including 0{t^), m'l is diagonal also up 
to and including C(t^), and the matrices W , W^, U and V are orthogonal up to the same 
order. These calculations have been checked by using two independent methods. 



Appendix C Notations and one-loop formulae for R 
and 6afy^^ 



IJ 



In this part we fix the notations and we report the formulae for Rij and for 5a^^^^ which 
we have used in section 5. The main references are [29-32] and we follow these. 
The mass matrix of the charginos is given by: 



Wr 



H 



+ h.c. 



IL 



with 



Mr 



M2 \/2miY cos P 
\^miY sin P fi 

This matrix is diagonalized by 2 x 2 rotation matrices Ol and Or as: 



OijM,0[ = diag(M~-,M~- 



(114) 



(115) 



(116) 



where the diagonalizing matrices connect mass and interaction eigenstates in the following 
way: 

' \ _ ^ \ ( XlR 



X2L 



' IL / \ X2R / \ -'-'2/? 

and the mass eigenstates are written as Xa — Xal + Xar (^ = 1)2) with masses M~-. 
The neutralino mass matrix is given by: 



H- 



(117) 



-Cn,D-i Bl Wl H'^j, H^j, ) Mm 



( B, \ 



/70 



+ h.c. 



^^The different factors of t and u are recovered in the final result by replacing x with j . 



:ii8) 



35 



with 



Ml 






Ma 



—niz sin 9w cos /3 sin sin /3 

cos 9w cos /3 — cos 9w sin /3 



— mz sin 9w cos /? cos 9w cos /9 
y sin 9w sin — cos 9w sin /? 



-/X 



We can diagonalize M^v by a rotation matrix O^r: 

OnMnOI = diag (m^, M^o, M^o, M^) , 
where Ojv connects mass and interaction eigenstates in the following way: 



/ 

(119) 
(120) 



X2L 
V X4L j 



O 



N 



( B, \ 

wl 

V ^2°L ) 



(121) 



and the mass eigenstates are given by = + (^ = 1, 2, 3, 4) with masses M^o . 
The mass matrices for the charged sleptons and for sneutrinos are given by: 



with 



m 



eLL 



m 



eLR 



^2 



m: 



^2 



'■eRL "^eRR 

where ^(ej/)LL ^"^^ ''^Irr i'^ general hermitian matrices and jfigj^^ = {p^'Rif ■ We 
diagonalize the mass matrix TWg by a 6 x 6 rotation matrix as: 



(122) 



(123) 



where the mass eigenstates are: 



'X,i+3''i 



(124) 



(125) 



with masses m~ {X — 1, ... ,6). 

•■JC 

Analogously, the sneutrino mass matrix is diagonalized by: 



where the mass eigenstates are: 



with masses m?^ {X — 1, 2, 3). 



(126) 
(127) 



The normalized branching ratios, Rij, for the LFV transitions I 

BR{k Ij-f) _ 487r3a 



Zj7 are: 



R 



(128) 



36 



and the decay rates are given by: 

r(/, ^ 1,1^,17,) = ^^m^ , r{h ^ /,7) = (l^^r + \A^n ■ (129) 

Each coefficient /Ig can be written as a sum of two terms: 

= + (130) 

where A^^^'^ and 742*^'^'^ stand for the contributions from the neutrahno and from the 
chargino loops, respectively. These coefficients are explicitly given by: 



ML ^ ^ L 

^ 327r2m? 



(131) 

and A^2^^ = A^^^Il^R with xax = /m? , and 

(132) 

and Af)'' = Af^^li^fi with = Myml^. 

The terms A'j^x and CiAx and the loop functions gin and gfjc read as follows: 

Nlkx = -4 \[-{0n)a,2 - {0^)A,itanew] u],^, + ^^(0^)^,3^ |(133) 

V 2 [ 771^^ COS p ) 

N^Ax = { ^^^^"^;^^ (Oiv)A3^i + 2(OAr)A,itan^^^C/t,^3| (134) 

= -y2(0i.)A,i^l,, (135) 

^^^x = 92 ^ A Ol)aM, (136) 

y Imw cos p 

and 



5'in(a;Ax) = 


6(1 


(1 6xax + ^x\x + '^^Ax Q^Ax In Xax) 

- XAxr 

X , ( 1 ^Ax + '^^Ax In Xax) 

- XAxr 


(137) 


g2n{xAx) = 


(1- 


(138) 


5'ic(a^Ax) = 


6(1 


(2 + 3xax Gx\x + + 62; AX In xax) 

- Xax) 


(139) 


g2c{xAx) = 


(1- 


( 3 + 4a; AX x\x 2 In a; ax) • 

XAxr 


(140) 



We note that the functions fin and fic , displayed in section 4, are related to the loop 
functions gin and gic, mostly through taking the first derivative. 

The deviation of the muon anomalous magnetic moment, 5a^^^^ , due to SUSY contri- 
butions can be written as: 

(n) (c) 
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with 



(n,c) (•n,c)L , (n,c)R 

9l -9l +9^ 



(142) 



These terms are exphcitly given by: 
1 mf 



Sc)L 



CiAxCiAx9lc{.X ax) + C^y^xCiAx9lc{xAx) + C^^X^iAX —92c{xax) 



and g'f'^ = gf^^\L^R with xax = M^./ml , and 



(143) 



Xa' 



ML 



1 

IGvr^ m? 



N!lxNtix9in{xAx) + Nj'j,xNjXxg^^{xAx) + NtAxNS,x^92n{xAx] 



(144) 

and = qI^^^Il^r with xax = M% /"^r • The terms NiAx and Cj^x and the functions 



gin and (^jc have been aheady introduced in eqs. (133 136) and in eqs. (137 140). 
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